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Abstract 

The phenomenon of forced unsteady separation 
and eruption of boundary-layer vorticity 1s a 
highly-complex, h1gh-Reynolds number flow phen- 
omenon, which abruptly leads to the formation of a 
dynamic stall vortex as demonstrated earlier by the 
authors for a NAGA 0015 airfoil undergoing con- 
stant rate p1tch-up motion. This, as well as the 
results of other researchers, have convincingly 
demonstrated a complex vortical structure within 
the state of unsteady separation prior to the evo- 
lution of dynamic stall. This phenomenon of vortex 
eruption, although observed In studying dynamic 
stall phenomena. Is also associated with transition 
from laminar to turbulence flow and Its generic 
nature has been stressed by many researchers 
Including the present Investigators. 

An unsteady Navler-Stokes (NS) analysis Is 
developed for arbitrarily maneuvering bodies using 
velocity-vorticity variables; this formulation Is 
nearly form-invariant under a generalized non- 
Inertlal coordinate transformation. A fully- 
ImpUdt uniformly second-order (except convective 
terms) accurate method Is used, with the nonlinear 
convective terms approximated using a biased third- 
order upwind differencing scheme to be able to 
simulate hIgher-Re flows. No explicit artificial 
dissipation Is added. The numerical method Is 
fully vectorized and currently achieves a computa- 
tional Index of 7 micro-seconds per time step per 
mesh point, using a single processor on the CRAY 
Y-MP 8/864 at the Ohio Supercomputer Center. The 
simulation results show that the energetic free 
shear from the leading edge Is responsible for the 
wall viscous layer to abruptly erupt near the cen- 
ter of the counterclockwise rotating eddy In the 
unsteady boundary layer. Primary, secondary, 
tertiary and quaternary vortices have been observed 
before the dynamic stall vortex evolves and gathers 
its maximum strength. This study will discuss the 
simulation results of Reynoldis number up to Re « 

45,000 and will also discuss the effort of Initial 
acceleration In a specific maneuver, on the evolu- 
tion of the stall vortex. 
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ABSTRACT 

Dynamic stall phenomena have been simulated numer- 
ically by a discrete vortex method and viscous flow calcula- 
tions. Both simulations could capture the characteristics of 
dynanuc stall phenomena qualitatively. Also some numeri- 
cal results show good agreements with experiments. 

In the calculations by a discrete vortex method com- 
bined with a panel method the potential flows around wing 
sections is expressed by vortex sheets and separated shezir 
layers are expressed by discrete vortices. Separated flows 
around pitching airfoils(NACA0012) are simulated for var- 
ious conditions. A hysteresis of lift of airfoil at dynamic 
stall is obtained. The results suggest that the method has 
a excellent capability of simulating vortical flows with ex- 
cellent small computation cost. In the calculations by vis- 
cous flow calculations incompressible Navier-Stokes equar 
tions have been solved by a third-order upwind scheme in 
order to understand the flow structure and mechanism of 
dynanuc stall. Especially the flow mechanics of movement 
of separation point on the wing surface, vortex formation 
from the surface of the wing and unsteady Kutta conditions 
at the trailing edge of the wing are investigated carefully. 
The unsteady flow fields around a pitching airfoil are cal- 
culated by moving a grid system relative to the freest ream. 
Remarkable chenacteristics of dynamic stall are obtained. 
Also the hysteresis curve of aerodynamic characteristics of 
Cl and Ci 4 are obtained. 

l^JNTRODUCTIQN 

Studies on dynamic stall have been conducted exper- 
imentaJly and theoretically by the mamy researchers*"'*^ 
However as the phenomena are essentially non-linear and 
the flow field becomes complicated as the separated region is 
increased, those phenomena are not understood sufficiently. 

In this paper two efforts for simulating dynamic stall 
phenomena are discussed. One is based on a discrete vor- 
tex method combined with a panel method and the other is 
based on viscous flow calculatioxis. Both simulations could 
capture the qualitative properties on dynamic stall phenom- 
ena. Also some numerical results show good agreements 
with experiments. 

In the calculations by a discrete vortex method com- 
bined with a panel method**^ *^ the potential flows around 
wing sections is expressed by vortex sheets and separated 
shear layers are expressed by discrete vortices. For the cal- 
culations a new numerical scheme developed by the present 
authors^ have been used. In the calculation a wing section 
is expressed by a set of linearly distributed vortex sheets and 


the separation points are determined by boundary layer cal- 
culations. Sepairated flows around pitching airfoils are sim- 
ulated for various conditions. A hysteresis of lift of airfoil 
at dynamic stall is obtained. The results suggest that the 
method has a excellent capability of simulating vortical flows 
with excellent small computation cost. 

In the calculations by viscous flow czJcuIations incom- 
pressible Navier-Stokes equations have been solved by a 
third-order upwind scheme in order to understand the flow 
structure and mechanism of dynamic stall. Especially the 
flow mechanics of movement of sepairation point on the wing 
surface, vortex formation from the surface of the wing and 
unsteady Kutta condition at the trailing edge of the wing 
are investigated carefully. The unsteady flow fields around 
a pitching airfoil are calculated by moving a grid system 
relative to the freestream. Remarkable characteristics of 
dynamic stall are obtained. Also the hysteresis curve of 
aerodynamic chajacteristics of C^, and Cm are obtained. 

2. Numerical Simulation of Dynamic Stall 
by Discrete Vortex Method 

2.1 Analytical Method 

The complex potential f of the flow field is expressed 
by the following form: 

/ = i^e-“2 + i jT -^log{z - Qd( 

( 1 ) 

, where U is the uniform flow velocity and a is the angle 
of attack. The discrete vortices expressing separated shear 
layers arc shed from the separation points on the wing sur- 
face, The circulation of each vortex is estimated using the 
velocity at the point assumed as the edge of the boundary 
layer near the separation point. Hence, unknown quantities 
in the flow field are the amount of strength of the singular 
points on the wing surface and solved at each time step by 
using boundary condition for norma] velocity through the 
control point on the surface. The detailed description of the 
method are discussed in Ref. 7, 

For the calculation of the flow around the pitching 
wing section, incident angle of the freestream is sinusoidally 
changed. Then the angle of attack of the wing section is 
given by following form: 



a = Ofo + aisinuft 


( 2 ) 


, where the angular frequency u is related to the reduced 
frequency k by the following form: 


k = 


UfC 


( 3 ) 


where c is the chord length. As the separation point on 
the wing section is not obvious a priori, the separation 
point is determined by the boundary layer calculations. The 
Thwaites* method* ' is used for the laminar boundary layer 
calculation and the Truckenbrodt^s method* » is used for 
the turbulent boundary layer calculation. 


2.2 Separated Flows around a Wing Section at a Fixed 
Angle of Attack^ 

Fig. 1 shows sevial calculated instantaneous stream 
lines around a wing section of NACA4412 at angle of at- 
tack of 20 degrees. As shown in the Figure the flow is quite 
unsteady and the changes of separation point and separated 
region can be observed. The calculated pressure distribution 
is comp«ired with the experiments in Fig. 2. As shown in 
the Figure the caJculated pressure distribution shows good 
agreement with that by experiments^ including the separ 
rated region. Fig. 3 shows Ci, and C© with a. A good 
agreement of Cl between the calculations and the experi- 
ments is obtained. Also the results of Co show good agree- 
ment with the experiments. Those results suggest that a 
discrete vortex method is quite powerful for prediction of 
aerodynamic characteristics of a wing section at static angle 
of attack. 




Fig. 2 Comparison of pressure distribution with experi- 
ments(NACA4412, a = 20degrees, i?e = 1.0 x I0«) 




Fig. 1 Instantaneous stream lines(NACA4412, a = 
20degreeSy Re ^l.Qx 10*) 


2.3 Separated Flows around an Oscillating Wing Section 

in. Pitch 

The cdculated results of sepwated flows around pitch- 
ing airfoil(NACA0012) under the condition of a = 15° -h 
10°sina/t(^ = 0.15) and Reynolds number=1.0 x 10® is 
shown in Fig. 4. As shown in the figure separated region 
is small in pitching-up process and it becomes much larger 
in pitching-down process. Quite different characteristics of 
flow patterns between in pitching-up and pitching-down pro- 
cesses are obtained. The results show some features of dy- 
namic stall. The calculated aerodynamic characteristics are 
compared with experiments^^ as shown in Fig. 5. The calcu- 
lated Hysteresis curve of Cl show almost same tendency as 
the experiments. The results show fairly good agreements 
with experiments. 

Another circulation have been conducted under the 
condition of a = 15°-h 14°sina;f(l: =: 0.1) and Reynolds num- 
ber=1.0 X 10® as shown in Fig. 6. In the Figure separated 
region is small in pitching-up process and it becomes much 


larger in pitching-down process. Quite different character- 
istics of flow patterns between in pitching-up and pitching- 
down processes are also obtained. The calculated aerody- 
namic characteristics are compared with experiments^^ as 
shown in Fig. 7. The calculated Hysteresis curve of 
show almost same tendency as the experiments. The results 
show fairly good agreements with experiments qualitatively. 
The results show of oscillating airfoil shows higher max^ 
imum value compared with that of a stationary airfoil. 

Those calculated results show qualitatively good agree- 
ment with experiments and excellent capability of the 
method is proved. Also as the computation time for each cal- 
culation is quite small compared with those of other finite- 
difference methods, the method is quite useful for the first 
estimation of the aerodynamic characteristics of a new wing 
section. 

3. Numerical Simulation of Dynamic Stall 
by Viscous Flow CalcrJation 

3_.l Numerical Procedures 

Unsteady flows of viscous incompressible fluid flow 
around a airfoil is considered. The governing equations are 
equations of continuity and incompressible Navier-Stokes 
equations: 

div V = 0 ( 4 ) 

dV 1 

-^ + ( V . V) V = -grad p 4- (5) 

where Re is Reynolds number. Following MAC method 
(Marker and Cell method)^°^ the Poisson equation for the 
pressure p is derived by taking divergence of Eq(5): 

Ap = -div(V . V) V -j- R (6) 

where 


Although R in Eq.(7) is identically zero due to Eq.(l), it is 
retained here as a corrective term in order to prevent the 
accumulation of numerical errors. If V is given at certain 
time, then the Poisson equation (6) can be solved to get 
p, and then by substituting these value into Eq.(5), V at 
next time is calculated from Eq.(5). The detailed numerical 
procedures will be given in Reference 11. 

For mesh grid generation Steger and Sorenson’s 
method^^^ is used. The method is useful to generate a body- 
fitted mesh system for an arbitrary body with keeping 
most orthogonal coordinates system at the wall. The mesh 
system is generated as 0-grid type and quite fine mesh is 
prepared in the vicinity of the body surface. For the calcu- 
lation of separated flow around oscillating airfoil a general 
coordinate system with time-dependent variables is intro- 
duced. However in the computation using those coordinate 
systems the displacement of the grid point and distortion of 
the mesh system for each time interval should be considered 
carefully. In the present computation the method proposed 
by Nakamichi^^^ is used. 

In the velocity fields non-slip conditions are used for 
the wing surface and uniform flow conditions are used for 
outer boundary. In the pressure fields zero derivative normal 
to the wing surface is used and zero derivative parallel to the 
freest ream is used at exit boundary. 

3.2 Separated Flows around an Airfoil at Fixed 
Attack Angle 

Before the calculations the preliminary calculation has 
been conducted for separated flows around circular cylinder 
at low Reynolds number in order to check the validity of 
the numerical scheme by- the present authors*^^ The results 
show excellent agreements with experiments. The width of 
the separated region and flow patterns are predicted quite 
precisely. The results show the present numerical scheme is 
quite useful and reliable for incompressible viscous flow. 

The mesh system(85 x 61) generated for a wing sec- 
tion of NACA0012 is shown in Fig. 8. Quite fine mesh 
is prepared in the vicinity of the wing surface and quite 
wide range of outer boundary are considered. The separated 
flows around airfoil at various attack angle is calculated at 
Reynolds number of 300. The representative results of in- 



Fig. 8 Mesh system (NACA0012;IMAX=85, JMAX=61) 



(a) a = 14® 


(c) a = 20® 

Fig. 9 Comparison of Flow patteins(instantaneous stream 
lines) between calculations and experiments(Aluminum 
powder) (NACA0012; Re=300) 



(b) cr = 16® 





stantajieous stream lines are compared with experimental 
results obtained by the present authors as shown in Fig. 9. 
The results show excellent agreements with experiments at 
various attack angle. As shown in the Figures separated 
region become large as attack angle is increased. In the 
present conditions separated shear layer forms a vortex and 
convects in the freestream. Karman vortex shedding pattern 
is observed behind an airfoil section. A separated shear layer 
shed from the trailing edge forms another vortex zmd rolls up 
behind a airfoil and convects downstream. The flow is quite 
periodic and those flow patterns repeat at every interval. 

The separated flow patterns can be also visualized by 
the electrolytic precipitation method. Fig. 10 shows com- 
parison of calculated electrolytic precipitation picture and 
that obtained by the experiments. In the calculation nu- 
merous marker particles aie introduced in the vicinity of 
the airfoil surface at some time interval and their temporal 
traces, which are induced by the velocity fields, are calcu- 
lated. Those traces form streak lines of the flow fields. A 
Crank- Nicolson^s implicit time integration is used in order 
to get accurate temporal resolution. The results show the 
numeric^ streak lines , i.e. numerical electrolytic precipi- 
tation pictures, is quite useful in order to understand the 
structure of the flow fields and compare the calculated re- 
sults with experiments. 


3.3 Separated Flows around an Oscillating Airfoil 

KI M ^ 

In the previous section the numerical scheme used in 
the present calculation is proved to be quite useful and reli- 
able, the unsteady flow fields around a oscillating airfoil in 
pitch are calculated by moving a grid system relative to the 
freestream. The flow conditions are selected from the ex- 
periments conducted by the present authors. The center of 
rotation is located at 25 percent of the chord and an airfoil 
section of NACA0012 is oscillated sinusoidally in pitch at 
reduced frequency o{ k = ufc/2U = 0.2 and Reynolds num- 
ber of 7.0 X 10 ^ c is a chord length and U is freestream 
velocity. 

The representative calculated results under the condi- 
tion of q; = 120-1 - k = 0.2 and Reynolds number=7.0 
X 10^ is shown in Fig. 12. Instantaneous stream lines and 
numerical smoke wire diagrams are shown in the Figures. 
For the calculations of numerical smoke- wire picture nu- 
merous fluid particles are introduced at equi-spacing points 
along a perpendicular axis, which is set in front of an airfoil 
for several time interval and their temporal traces, which are 
induced by the velocity fields, are calculated. Those traces 
form streak lines and a numerical smoke- wire picture is ob- 
tained. A Crank-Nicolson^s implicit time integration is used 



time = 84.82 a = 15.80* up 


time = 87.96 a = 15.80* down 





time = 81.68 a = 14.35* up 


time = 91.11 a = 14.35* down 



time = 78.54 a = 12.00* up 


(a) Instantaneous stream lines 

Fig. 11 Calculated flow patterns(NACA012, a =129 + 4<>sinwt(Jt = 0.2), Re = 7.0 x 10^ 





Fig. 11 (b) Numerical smoke- wire picture 


in order to get accurate temporal resolution. The results 
show the numerical smoke-wire pictures are quite useful in 
order to understand the structure of the flow fields and com- 
pare the calculated results with experiments. 

As shown in the figure separated region is small in 
pitching-up process and it becomes large in pitching-down 
process. Quite different characteristics of flow patterns be- 
tween in. pitching-up and pitching-down processes are ob- 
tained. The significant change of the flow patterns in 
pitching- up and p itching-down process show the e xistence of 
the Hysteresis of the aierodynamic characteristics. The hys- 
teresis curve of aerodynamic characteristics of Cl and Cm 
at the same condition are shown in Fig. 13. In particular 
Cl at higher attack angle than the stalling angle of this air- 
foil section keeps higher value, then decreases rapidly. The 

results shows remarkable characteristics of dynamic stall. 

4. CONCLUSIONS 

Dynamic stall phenomena have been simulated numer- 
ically by a discrete vortex method and viscous flow calcular 
tions and major results of using both numerical efforts are 
discussed. The major conclusions of the present study are 
summarized as follows: 

In the calculations by a discrete vortex method com- 
bined with a panel method the potential flows around wing 
sections aerodynamic characteristics of dynamic stall is cal- 
culated properly and a hysteresis of lift of airfoil at dynamic 
stall is obtained. The results suggest that the method has a 
excellent capability of simulating vortical flows with excel- 
lent small computation cost. 


In the calculations by viscous flow calculations sepa- 
rated flows around oscillating airfoil in pitch are simulated 
by using the third-order upwind scheme and a moving mesh 
system. The calculated separated region is small in pitching- 
up process auid it becomes large in a pitching-down process. 
Quite different characteristics of flow patterns between in 
a pitching-up and pitching-down processes and hysteresis 
curve of aerodynamic characteristics are obtained. 
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Fig. 12 The hysteresis curves of aerodynamic characteristics of Ci, and C*, (NACA012 a 
12® 4- 4»sin£j<(i = 0.2),Re = 7.0 X 10<) 
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Fig. 7 Comparisons of calculated hysteresis curves of Ci with experiments^ ( o = 15 
14»s»nw<(it = 0.1; NACA0012, Re = 1.0 x 10*)) 
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